Anomalous pinning behavior 
in an incommensurate two-chain model of friction 
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Pinning phenomena in an incommensurate two-chain model of friction are studied numerically. The 
pinning effect due to the breaking of analyticity exists in the present model. The pinning behavior is, 
however, quite different from that for the breaking of the analyticity state of the Frenkel-Kontorova 
model. When the elasticity of chains or the strength of interchain interaction is changed, pinning 
force and maximum static frictional force show anomalously complicated behavior accompanied by 
a successive phase transition and they vanish completely under certain conditions. 



I. INTRODUCTION 

In recent years, the study of friction has been attract- 
ing much attention in physics Q|. Nano-scale frictional 
phenomena have been examined experimentally using 
frictional force microscopes quartz micro- valance 
techniques || and so on. In theoretical studies, the 
Frenkel-Kontorova (FK) model Q| and its related ones 
have been employed as a promising model of such nano- 
scale friction by several researchers |3|-|| . The FK model, 
in general, consists of an atomic chain on a substrate with 
periodic potential. In the chain harmonic force works 
between neighboring atoms. When the mean atomic 
distance and the period of the potential is incommen- 
surate, the FK model shows a phase transition, which 
has been discussed in detail by Aubry and coworkers 
P,^0|. Hence this phase transition is called the Aubry 
transition. The Aubry transition has following features. 
When the amplitude of the substrate potential is smaller 
than a certain critical value, the lowest phonon exci- 
tation is gapless and therefore a free sliding mode ap- 
pears. This means vanishing maximum static frictional 
force. Above the critical amplitude, however, the atoms 
in the chain are pinned strongly nearby the potential 
minima and a finite gap exists in the phonon excita- 
tion. This state is called the breaking of the analyticity 
state. Then finite energy is needed to slide the chain, and 
therefore the maximum static frictional force becomes fi- 
nite. The extended FK model which consists of interact- 
ing two deformable chains also have been investigated so 
far. The static structural properties of two-chain models 
have been investigated in Refs. 11 where each chain 
is often treated as a continuum elastic line. The con- 
tinuum approximation works effectively in the study on 
the commensurate- incommensurate transition | fl2| . How- 
ever, the pinning effect, that arises from the discrete 
nature of lattices, are smeared out inevitably. In other 
words, the Aubry transition never occurs in the contin- 
uum models. On the basis of a two-chain model with 



discrete lattice structures, Matsukawa and Fukuyama 
investigated the static and the kinetic frictional forces 
p5|,^6[ . The model proposed in their study consists of 
two atomic chains, where interchain atomic force works 
between atoms in one chain and in another and harmonic 
force works between neighbor atoms in each chain. In 
some cases of elastic parameters of chains, they found 
that the maximum static frictional force for the two-chain 
model becomes larger than that for the FK model with 
the same strength of the interchain interaction. Fur- 
thermore, they discussed the relationship between the 
strength of the maximum static frictional force and the 
velocity dependence of kinetic frictional force. 

In this paper, we revisit the two-chain model of fric- 
tion employed in Ref. JlJ] and examine the frictional phe- 
nomena in a wide range of model parameters. In partic- 
ular, pinned states are investigated thoroughly in con- 
nection with the breaking of analyticity state due to the 
Aubry transition. It turns out that the maximum static 
frictional force shows complicated behavior against the 
change in elastic parameters and vanishes completely in 
certain conditions. This anomalous pinning behavior is 
discussed in relation to the static lattice structures. We 
also focus on the velocity dependence of the kinetic fric- 
tional force in sliding states. 



II. TWO-CHAIN MODEL OF FRICTION 

The two-chain model of friction employed here is sum- 
marized in the following |l5|j . We consider two atomic 
chains, i.e., an upper chain and a lower chain. Each 
atom has one-dimensional degree of freedom parallel to 
the chain. Intrachain interaction with harmonic form 
and interchain interaction are taken into consideration. 
The effects of energy dissipation are assumed to be pro- 
portional to the difference between the velocity of each 
atom and that of the center of gravity of the chain. The 
upper chain is driven by the external force parallel to the 



1 



chain. Assuming overdamped motion, we get the equa- 
tions of motion of the atoms in the upper and the lower 
chains given by 

m a j a (ui - (v,i)i) = K a (u i+ i + Uj_i - 2m) 

+ ^F I {u i -v j )+F ex , (1) 
jeb 

m b j b (vi - (vi)i) = K b (v. l+ i + Vi-x - 2vi) 

+ Y^Fi{v i -u j )-K s {v i -ic b ), (2) 

j£a 

where m («,-), m a (mi), 7 a (7b), K a (Kb), and N a (N b ) 
are the position of the ith atom, the atomic mass, the 
parameter of energy dissipation, the strength of the in- 
teratomic force, and the number of atoms in the upper 
(lower) chain, respectively. K s denotes the strength of 
the interatomic force between the lower chain and the 
substrate, which is necessary to bind the lower chain. ( ) ; 
represents the average with respect to i. Fj and F ex are 
the interchain force between the two atomic chains and 
the external force, respectively. The interchain atomic 
potential is chosen as follows: 



Ui = — — exp 



(3) 



where Kj is the strength of the interchain potential, and 
c b the mean atomic spacing of the lower chain. The in- 
teratomic force is given by Fi(x) = —4-Ui. The time- 
averaged total frictional force of the present model is 
given by 
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fric 
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(4) 
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This expression of the frictional force is valid for both 
static and kinetic ones. The form of interchain force, 
Fj(x) = —4-Ui, is approximated in the following two 
limiting cases JTtJ . When the atoms in the lower chain are 
rigid and fixed at the regular sites, which corresponds to 
the case of Kb/K a — > oo or K s /K a — ► oo, the interchain 
force that acts on the upper chain is approximated by 
one term in a Fourier series: 



jeb 



-0.47-ffjsin 



Ui 

2tt— 

Cb 



(•5) 



Vj=JC b 

Then the two-chain model is reduced to the FK model. 
In the opposite limit, when the upper chain is fixed at 
the regular sites (K a /K b 1), the interchain interaction 
that acts on the lower chain is also described by one term 
in a Fourier series: 



n„ 



^Fiiui-Vj) 



-0.83K T sin 2tt 



(6) 



where c a is the mean atomic spacing of the upper chain. 
These approximations on the interchain force are valid 
only in the above two limiting cases on elastic parame- 
ters, and they may break in intermediate cases. This is 
a crucial point in the following discussion. 



III. NUMERICAL METHOD 

For the numerical simulation of the model, the Runge- 
Kutta (RK) formula is employed to solve the equations of 
motion. The periodic boundary conditions are employed 
in both chains. Hence, the ratio c a /c b is equal to N b /N a , 
where N a and Nb are the numbers of atoms in the upper 
and lower bodies. 



- = — = — = 1.618. 
144 



Ca 
Cb 



N b 
N 



(7) 



where the ratio is determined by using a continued- 
fraction expansion of the golden mean to emulate incom- 
mensurability. Throughout the present study we set the 
values of the model parameters as 



N a = 144, N b = 233, c a = 1.618 • 

K a = 1,7a = 76 = 1- 



•,c b = l,m a = m b = l, 
(8) 



We control the elasticity of chains by changing the spring 
constants K s and K b , and the strength of interchain in- 
teraction by Kj. 

In numerical simulations we mainly employ an initial 
atomic configuration where atoms are located at regular 
sites periodically. During the RK steps the chains are 
relaxed to stable states and finally reach there. In Section 
IV-A-5 we will refer to the effect of the initial atomic 
configuration in connection with pinning behavior. 

We employ the following numerical criterion and meth- 
ods in the calculations of several quantities. If iii and Vi 
calculated during the RK steps satisfy a velocity con- 

dition, v / (E^0 2 + £fW)/W 1 + A r 6 ) < io- 10 , 

the RK calculation is stopped, and the state obtained 
then is considered to be static. The phonon frequency is 
calculated using a dynamical matrix for this stationary 
state. The maximum static frictional force is evaluated 
as the critical force above which the velocity condition 
is not satisfied. Using the criterion, the difference be- 
tween pinned and sliding states is distinguishable. The 
frictional force is calculated according to Eq. (4) or the 
method used in Ref. |l5|. In a sliding state, after the 
system reaches a steady state, the temporal average is 
performed in calculating the kinetic frictional force over 
a time period much longer than a time during which 
the center of gravity of the system moves by the system 
length, N a c a (= N b c b ). 
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IV. RESULTS 
A. Pinned states 

1. Lowest phonon frequency 

In this section we investigate the pinning effect of the 
two-chain model in the absence of external force. To 
investigate the feature of pinned states, we first calcu- 
late the lowest-phonon frequency, which is a significant 
quantity because finite lowest-phonon frequency, i.e., the 
phonon gap, means the presence of pinning and its square 
is proportional to the restoring force due to pinning ef- 
fects. 

For the comparison with the two-chain model, we first 
show the squared lowest phonon frequency {ujf pi ) of the 
FK model described by Eqs. (1) and (5) with Kj = 1 in 
Fig. 1 . In this case the elastic parameter is K a only, and 
there exists its critical value. Below the critical K a , uif p{ 
becomes finite, corresponding to the appearance of the 
breaking of analyticity state due to the Aubry transition. 
Above the critical K a , uf pi vanishes. The change ofu)f pi is 
continuous at the critical point. Such behavior of ujf p{ for 
the FK model was reported in Ref. JlOj] . In the present 
study on the two-chain model we have another elastic 
parameters in the lower chain. Hence, we focus on the 
effect of the elastic relaxation of the lower chain, and then 
K s and Kb are chosen as variable elastic parameters. 

In Fig. 2(a) we show uif p{ as a function of an elastic 
constant Kb in the case of strong interchain interaction 
Ki = 1, and K s = 1. The strength of the interchain 
interaction (Kj = 1) is chosen to be greater than the 
critical value of the Aubry transition for the FK model, 
KYyy^ 1 ~ 0.33 for K a — 1. Because of strong interchain 
interaction, the phonon gap is finite and almost constant 
in the region of large value of Kb (2.4 < Kb). However, 
steep valleys appear in the range 0.81 < Kb < 2.4. The 
amount of u>f { changes there more than two orders of 
magnitude. It is considered that utf pi vanishes at each 
valley. The finite values of ujf p{ at the bottoms of the val- 
leys will be due to the numerical accuracy of the present 
calculation and the finite magnitude of changing in Kb- 
When Kb becomes smaller than a certain value (~ 0.81), 
ujf pi increases sharply and then becomes almost constant. 
The vanishing phonon gap at each valley seems to indi- 
cate a sort of phase transition. 

When the interchain interaction Ki is somewhat weak- 
ened, but its strength is still greater than K"p^ al , more 
drastic behavior of phonon gaps is observed. Figure 
2(b) shows uf pi obtained for Kj = 0.45 and K s = 1. 
Finite phonon gaps exist both in small and large Kb 
regimes, and steep valleys appear in the intermediate 
regime, 0.2 < Kb < 0.6. Such behavior is quite simi- 
lar to that for Kj — 1. In a wide regime, 0.6 < Kb < 1.4; 
however, the phonon gap vanishes completely within a 



numerical accuracy. This regime is quite distinct from 
other regimes with narrow valleys seen in Fig. 2(a) as 
discussed in the next subsection. The anomalous behav- 
ior of the phonon gap disappears when a large value of 
K s is chosen as in Fig. 2(c), where K s — 10 (Ki = 1). It 
is obvious that the behavior of the phonon gap depends 
on the elastic parameter K s as well as Kb- In Fig. 2(a) 
the behavior of the phonon gap looks self-similar. If this 
is the case, the phonon gap reveals more complicated be- 
havior against a smaller change in Kb in the intermediate 
Kb regime. 

The discontinuous and complicated behavior of ujf p{ 
against the change in elastic parameters for the two-chain 
model is obviously quite different from that for the FK 
model shown in Fig. 1. We confirmed that the magnitude 
of the phonon gap at all KbS chosen in Figs. 2(a)-(c) is 
entirely insensitive to the enlargement of the system size 
determined by using a continued-fraction expansion of 
the golden mean in Eq. (7). 

2. Maximum static frictional force 

Next we calculate the maximum static frictional force 
numerically by applying the external force to the pinned 
states. Figures 3 (a) and (b) show the maximum static 
frictional force calculated for pinned states shown in Figs. 
2(a) and (b), respectively. The maximum static frictional 
force also shows anomalous behavior, which obviously re- 
flects the phonon gap structures in Figs. 2(a) and (b). 
In Fig. 3 (a), where the values of the parameters, Kj = 1 
and K s = 1, are the same with those in Fig. 2(a), the 
maximum static frictional force shows multivalley struc- 
tures. The magnitude of the maximum static frictional 
force is finite in the whole Kb regime. We note that the 
values of Kb at the local minima and maxima of the max- 
imum static frictional force do not correspond exactly to 
those of W[ 2 pf shown in Fig. 2(a). This may be considered 
to be the effect of the external force, by which the pinned 
lattice structures are distorted and the depinning thresh- 
old force would be affected slightly. On the other hand, 
as seen in Fig. 3(b) for Kj — 0.45 and K s = 1, the max- 
imum static frictional force vanishes completely in the 
characteristic Kb regime where the completely vanishing 
phonon gap is observed in Fig. 2(b). In the case that 
Ki = 1 and K s = 10, it is confirmed that the maximum 
static frictional force as well as the phonon gap in Fig. 2 
(c) does not show any anomalous behavior and changes 
smoothly with Kb- 

3. Hull functions and lattice structures 

To investigate further the above anomalous pinning 
behavior observed for the phonon gap and the maximum 
static frictional force, we analyze the lattice structures of 
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the pinned states by examining hull functions. Although 
a hull function has been employed to analyze the breaking 
of analyticity state for the FK model [6 10 1, it has been 
reported that hull functions defined in two chains are also 
useful to analyze the lattice structures both in pinned 
and sliding states for the two-chain model |L7| . The hull 
functions for the two-chain model arc defined as 



Ui — ic a + a + h a (ic a + a), 
Vi = ic b + + h b (ic b + 0), 



(9) 
(10) 



where h a and hb are the hull functions in the upper 
and lower chains, respectively, and a and f3 are constant 
phases. The periodicities of the hull functions are ex- 
pressed as 

h a (x) = h a (x + Cb), h b {x) = h b {x + c a ). (11) 

When the chains are not deformed and hence the atoms 
are arrayed periodically, the interchain interaction po- 
tential is sinusoidal as mentioned in Eqs.(5) and (6), and 
then the position of the potential maximum in one pe- 
riod is located at the half of the period of the hull func- 
tion, x = Cb/2 (c a /2) for h a (x) (hb(x)) in our choice, 
a = P = 0. 

For the convenience of later discussions, we briefly 
summarize here some features of the hull function for the 
FK model. When the strength of the interchain interac- 
tion is less than the critical value of the Aubry transition, 
the hull function is smooth and continuous. Above the 
critical point, however, the breaking of analyticity due to 
the Aubry transition occurs, and then the hull function 
changes its form from continuous to discrete and shows a 
complicated structure with many gaps. Among the gaps 
the largest one is located at the half of the period of the 
hull function. The continuous form means continuous 
spatial atomic distribution in the underlying potential, 
and then no gap exists in the phonon excitation. On 
the other hand, the discrete one corresponds to a pinned 
state, which is accompanied by a finite gap in the phonon 
excitation. The spatial atomic distribution is vanishing 
at the maxima of the potential and the atoms are con- 
fined nearby the minima of the potential. Then the hull 
function shows the largest central gap, which character- 
izes the breaking of analyticity state for the Aubry tran- 
sition. 

Now we consider the case of the two-chain model. Fig- 
ures 4 (I)- (VII) show the hull functions h a and hb for sev- 
eral values of Kb indicated by arrows (I)- (VII) in Fig. 
2(a), where the magnitudes of the parameters, Kj = 1 
and K s — 1, are the same as those in Fig. 2(a). The 
hull function h b in Fig. 4(1) (K b = 0.1; (I) in Fig. 2(a)) 
shows the largest central gap at x = c a /2 « 0.809. This 
gap structure is essentially the same with that for the 
Aubry transition in the FK model and indicates that the 
lower chain is in the conventional breaking of analyticity 



state due to the Aubry transition. On the other hand, 
h a shows a discrete feature but does not have a central 
gap. Such a state in the upper chain where the central 
gap of the hull function is absent is not well defined in 
the context of the conventional breaking of analyticity 
state due to the Aubry transition of the FK model, but 
it is obviously a sort of breaking of analyticity states be- 
cause of the presence of gaps of the hull function. These 
gap structures of h a and h b remain even for Kb = 0.702 
(Fig. 4(H)). This reflects the constant ujf pf observed in 
the small K b regime (K b < 0.81) in Fig. 2(a). As K b 
increases further [K b > 0.81), however, the central gap 
of hb is destroyed and no central gap exists both in h a 
and hb- Some other gaps also vanish or shrink, otherwise 
enlarge, and furthermore, new gaps appear at several po- 
sitions. Figures 4(III)-(V) show h a and h b at several val- 
ues of Kb where ujf pi shows a local maximum against the 
change in Kb (see arrows (III)- (V) in Fig. 2(a)). It should 
be noted here that the gap structures of h a and hb shown 
in (III)-(V) are different from each other. In a Kb regime 
between two nearest- neighbor valleys of uif { , both the 
gap structures of h a and h b are almost unchanged. Only 
when K b is changed crossing through the valley of W[ 2 pf , 
the gap structures are suddenly changed, i.e., new gaps 
appear. In the large K b regime (K b > 2.4), as shown 
in Fig. 4(VI), the central gap appears in h a , but it is 
absent in h b , and the gaps of h b as a whole are highly 
reduced. This behavior indicates that the upper chain is 
in the conventional breaking of analyticity state due to 
the Aubry transition. These gap structures of h a and h b 
retain up to infinite K b while the amplitude of h b as a 
whole shrinks as K b is increased (Figs. 4(VI) and (VII)). 
The ujf pi also does not change in this large K b regime. 
To see how the change of gap structures of hull functions 
takes place at the valleys of uS^ pi , in Figs. 5 (a) and (c), 
we show hull functions for two nearest-neighbor valleys 
(local minima) of ojf pi at K b = 1.18 and 1.58. Note here 
that the K b for Fig. 4(111) is located in the K b regime 
between these two nearest- neighbor valleys of ojf pi . For 
the comparison with these states at valleys, the graph 
of Fig. 4(111) is plotted again in Fig. 5 (b). In Figs. 5 
(a) and (c) several gaps of the hull functions observed 
in Fig. 5(b) are destroyed by the appearance of certain 
states in gaps, i.e., the formation of new gap structures at 
the Kb's. As mentioned above the gap structures of hull 
functions for Fig. 5(b) are stable and almost unchanged 
in the regime of 1.18 < Kb < 1.58 where no valley of u>f { 
exists. Similar behavior is observed around each valley of 
cj[ 2 pf . When Kb reaches one of the critical value, the old 
gap structure becomes unstable and new states appear in 
gaps, which accompany the decrease of the phonon gap 
and then the pinning force. When Kb crosses the critical 
value, new gap structure becomes stable and the phonon 
gap increases. Further change of Kb moves the system 
to the next critical value and then the successive phase 
transition occurs. 
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In Fig. 6 we show the hull functions h a and hb for 
several values of Kb in the case of weak interchain in- 
teraction, K] = 0.45, and K s = 1, which are the same 
values with those in Fig. 2(b). In the small and large Kb 
regimes ((I), (II), and (VI), (VII) in Figs. 6, respectively) 
where uif { is almost constant, both the gap structures of 
h a and hb are unchanged, one of which shows the cen- 
tral gap. This indicates that one of the upper and lower 
chains is in the conventional breaking of analyticity state 
due to the Aubry transition in the large and small Kb 
regime, respectively. In the intermediate regime (Figs. 
6(III)-(V)), the gap structures as a whole are rather sen- 
sitive to the change in K b , but also in this case the gap 
structures are almost unchanged in a small Kb regime 
between two nearest-neighbor valleys of uif pi . Thus, the 
behavior of h a and hb is similar to that in Fig. 4 in 
this regime. However, in the Kb regime 0.6 < Kb < 1.4, 
for Kj = 0.45, where ojf { is vanishing, both the hull 
functions show a peculiar feature. In Fig. 7 we show 
several typical h a and hb in this Kb regime. It is clearly 
observed that both hull functions h a and hb are contin- 
uous and show sinusoidal forms, which are quite similar 
to that of the hull function for the FK model in the ab- 
sence of the breaking of analyticity. These correspond to 
states in which all atoms of both chains locate near its 
regular sites periodically and are weakly affected by the 
almost sinusoidal interchain force caused by atoms in the 
other chain. Since the continuous hull functions mean 
that the atomic distribution is spatially continuous both 
in the upper and lower chains, every atom in the upper 
chain moves smoothly when the upper chain is driven by 
the external force. Therefore, there are no energy costs 
against the sliding motion of the upper chain. Hence, the 
maximum static frictional force vanishes as observed in 
Fig. 2(b). 

For K s = 10 and Ki = 1 (Figs. 8(1)- (IV)), K does not 
show any remarkable changes of gap structures, hb grad- 
ually changes its gap structure only in a large Kb regime 
(Kb > 10), but the effect of the change is almost negli- 
gible because the amplitude of the gaps of hb becomes 
very small for such large K^s. Therefore all elastic ef- 
fects come from the upper chain. The behavior of the 
hull functions reflects the smooth change in the phonon 
gap in Fig. 1(c). In all regime of Kb, h a shows the largest 
central gap, but hb does not. Thus, in these pinned states 
the conventional breaking of analyticity due to the Aubry 
transition occurs in the upper chain for the whole range 
of Kb- Then the magnitude of the central gap in h a is 
quite larger than those of the gaps in hb and almost un- 
changed against the change in Kb- 
it can be confirmed that the breaking of analyticity 
states exist in the present two-chain model, but they are 
rather complicated and different from the conventional 
one due to the Aubry transition of the FK model. In the 
small (large) Kb regime, however, the gap structure of hb 
(h a ) is the same as the conventional one observed in the 



breaking of analyticity state in the FK model, while that 
of h a (hb) is different. That is, in the two-chain model, it 
is considered that whether the conventional breaking of 
analyticity state characterized by the largest central gap 
exists in the upper or lower chain depends on the elas- 
ticity of the two chains. When the lower (upper) chain 
is highly stiffer than the upper (lower) one, i.e., Kb/K a 
or K s /K a ^> 1 (<C 1), the atoms in the upper (lower) 
chain tend to relax into potential minima created by the 
atoms in the lower (upper) chain, and then the conven- 
tional breaking of analyticity state appears in the upper 
(lower) chain. In the intermediate Kb regime, however, 
quite different states from the conventional breaking of 
analyticity state observed for the FK model appear. We 
will again discuss this point later by calculating the en- 
ergy quantities of the system. 

It should be noted here that all the hull functions 
shown in Figs. 4-8 do not contain irregular points that 
break a rotational symmetry by tt of hull functions. This 
fact means that atomic configurations obtained above are 
not disordered, but they exactly reflect the discreteness 
of hull functions in an incommensurate system. This fea- 
ture of the pinned atomic configuration is the same as 
that for the FK model @. 

It is helpful here to observe the change of the pinned 
lattice structures in real space. Figure 9 shows the lo- 
cal lattice structure in the pinned state for Kj — 1 and 
K s — 1, which are the same values with those in Figs. 
2(a), 3(a) and 4. Here the atomic displacements from 
the regular periodic sites in the two chains Sui and Svi 
are plotted in Figs. 9(a) and (b), respectively for Kb's 
indicated by arrows (II)-(VI) in Fig. 2(a). In the small 
Kb regime (Kb < 0.8), the lattice structure is essentially 
unchanged (Figs. 9(a)-(l) and (b)-(l)). In the interme- 
diate regime (0.8 < Kb < 2.4), however, both Sui and Svi 
are very sensitive to the change in Kb- Figures (2) and 
(3) in both of Figs. 9(a) and (b) correspond to atomic 
displacements at local maxima of the phonon gap ojf pi 
in Fig. 2(a). It is obvious that the spatial modulation 
patterns of Sui and 5vi show quasiperiodicity approxi- 
mately for each Kb, but the spatial patterns are different 
for different values of Kb- When Kb increases further 
(Kb > 2.5), the reconstruction does not occur any more 
and the lattice structures as observed in Figs. 9(a)-(5) 
and (b)-(5) retain up to infinite Kb, but the atomic dis- 
placement in the lower chain Svi as a whole decreases its 
magnitude more and more. Similar changes of local lat- 
tice structures are observed also in the intermediate Kb 
regime for K] — 0.45 and K s = 1. Note here that if the 
behavior of the phonon gap has a self-similar nature as 
noticed in Fig. 2(a), then infinite sorts of local lattice 
structures would exist in the intermediate Kb regime. 
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4- Analysis of energy and discussion on the pinning 
mechanism 

We discuss further the anomalous pinning behavior by 
examining the energy of the system. Here it is useful to 
define the following energy quantities. Elastic energy in 
the upper chain -E e ia-a and that in the lower chain .Eeia-b 
are given by 

K a Na 

£cla-a = ~Y ^2(u l+1 - Ui - C a ) 2 , (12) 

K N " 

£cla-b = ~Vi -C b ) 2 . (13) 

i 

Elastic energy between the lower chain and the substrate 
£ i a _ s is given by, 

K Nb 

£da-s = -y £>i - lC b ) 2 . (14) 

i 

Total elastic energy i^ia-totai is the sum of 

Sola-total = -Eola-a + ^ela-b + ^cla-s- (15) 

Interchain interaction energy E- mt is expressed as 

1 N a N h 

E ™t = 2 \ Ul i u i ~ v i) ~ U i( ic a - 3Cb)\ , (16) 

i 3 

where the contribution in the case of periodic rigid atomic 
configurations is subtracted. To evaluate Eqs.(12)-(16) 
we used the atomic configuration obtained in the cal- 
culation of a;^. Figures 10 (a)-(d) show these energy 
quantities plotted against Kb for Kj = 1 and K s = 1, 
which are the same values with those in Figs. 2(a), 3(a), 
4, 5 and 9. When Kb is very small (Kb « 1), the ab- 
solute value of the interchain interaction energy |2?i n t| is 
much greater than i? e i a _ a , but not so much greater than 
-E^eia-s- The upper chain deforms little, but the lower 
chain deforms so large to gain the interchain interaction. 
In fact the hull function of the lower chain hb shows the 
central largest gap, but h a does not and its magnitude 
is small. These behaviors indicate that the upper chain 
has an almost periodic structure, but the lower chain ad- 
justs its local period to that of the upper chain to gain 
the interchain interaction, that is, the lower chain forms 
a kind of discommensurate structure or the soliton lat- 
tice. When Kb is increased, the deformation of the lower 
chain and -E e ia-s decrease. On the other hand, -Eeia-b 
increases because its coupling constant, itself, increases. 
When i^int becomes comparable to the total elastic en- 
ergy, the lattice structures of both chains in the small 
Kb regime becomes unstable and another structure ap- 
pears. The presence of kinks for all energy curves in 
Fig. 10 indicates that the structural change is a phase 



transition of the first order. The difference in the lattice 
structures between phases is observable directly in Figs. 
9(a) and (b). The disappearance of the central gap of the 
hull function h b in Fig. 4 is attributed to the structural 
phase transition. Both chains form complicated discom- 
mensurate structures where both chains deform to make 
the local commensurate structure in order to gain the 
interchain interaction. The occurrence of the structural 
phase transition leads to the sudden decrease (or vanish- 
ing) of the phonon gap at Kb « 0.81 as observed in Fig. 
2(a). In the intermediate Kb regime (0.81 < Kb < 2.4), 
many, or an infinite number of lattice structures appear 
there, each of which corresponds to a different discom- 
mensurate structure. When K b is increased small or 
infinitesimally, lattice structures change to another one 
by a structural phase transition. Such phase transitions 
occur successively in this regime against the change in 
Kb- As the lower chain becomes stiffer, the atomic dis- 
placement in the lower chain 5i>i is suppressed, and then 
-Eeia-b and .Eeia-s decrease. Because in this Kb regime a 
slight change in K b causes a structural phase transition, 
the phonon gap and the hull functions also change their 
structures correspondingly. It should be noted that the 
absence of the central gap in both hull functions char- 
acterizes the lattice structures in this intermediate Kb 
regime, where the successive structural phase transition 
occurs (Figs. 4(III)-(V)). This indicates complicated dis- 
commensurate structures of both chains. In the large Kb 
regime (Kb > 2.4), the structural phase transition docs 
not occur and a quite stable lattice structure appears 
in each chain, which is characterized by the largest cen- 
tral gap of h a . These behaviors indicate that the lower 
chain forms an almost periodic structure, but the up- 
per chain forms a discommensurate structure. In this 
regime the major contribution to -E e ia-totai comes from 
-Ecia-a because the atomic displacement 5vi is suppressed 
and SuiS are fixed into a locally commensurate configu- 
ration. Therefore, E e \ a _ a and E- 1tA are almost constant. 
Consequently, almost the same gap structures of the hull 
functions as those in Fig. 4(IV) and the same atomic 
configurations as those in (5) in Figs. 9(a) and 9(b) hold 
up to infinite K b . We note here that these features, such 
as the kinks for energy and the vanishing of the phonon 
gap, which mean a first-order phase transition, are never 
expected in the case of the FK model because the Aubry 
transition observed in the FK model is a higher-order 
phase transition. 

In Fig. 11 we show the energy quantities for the pa- 
rameters Kj = 0.45 and K s = 1 employed in Figs. 2(b), 
3(b), 6, and 7. Also, in this case kinks are observed for all 
energy curves. Within the Kb regime (0.6 < Kb < 1.4) 
where the phonon gap and maximum static frictional 
force are vanishing completely, any anomaly such as a 
kink is not observed. It is considered that one particular 
lattice structure of both chains appears in this K b regime 

The energy quantities obtained for K s = 10 and Kj = 
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1 are plotted in Figs. 12 (a)-(d). In this case, because 
of the large value of K s , the lower chain is much stiffer 
than the upper chain. Then all energy quantities change 
smoothly with Kb- This means that these pinned states 
in all Kb regimes are understood essentially in the context 
of the conventional breaking of analyticity state due to 
the Aubry transition for the FK model. Namely, the 
feature of pinned states for K s = 10 and K] = 1 are the 
same as those in the large Kb regime shown in Fig. 10, 
where K s = 1 and Ki = 1. 

For simplicity, we have changed the elasticity of the 
lower chain, while the elasticity of the upper chain is 
fixed (i.e., K a = 1). In general, it is considered that 
the relevant parameters to the pinning behavior of the 
present model are normalized ones: Kb/K a , K s /K a , and 
Ki/K a . 

5. Hysteretic behavior of phonon gap 

We have discussed the structural phase transitions of 
the first order. Then we expect hysteretic behaviors 
around the critical points. It is to be noted here that 
in all the results shown in Figs. 1-12, the initial config- 
urations of atoms for each value of Kb are regular and 
periodic. In Fig. 13 we show phonon gaps for the same 
parameters (K] = K s = 1) as those in Fig. 2 (a) when 
Kb is swept, that is, the initial atomic configuration for a 
certain value of Kb is the final stable configuration for its 
previous value. The solid lines with circles indicate the 
results obtained for the sweep with (a) increasing and 
(b) decreasing value of Kb- For comparison, the results 
in Fig. 2(a) are also plotted by the dashed lines with 
squares. By comparison, between these two results, the 
hysteretic behavior of the phonon gap is obviously ob- 
served in the intermediate Kb regime. Such hysteretic 
behavior, which appears suddenly only in the interme- 
diate Kb regime, is another evidence of the first-order 
phase transition. 

We have considered only three initial conditions corre- 
sponding to a periodic atomic configuration (Figs. 2-12) 
and two history-dependent configurations (up and down- 
ward sweeps of Kb) (Fig. 13). The pinning behavior 
for other initial conditions has not been investigated ex- 
haustively because a great variety of configurations may 
be considered. We have found here that the hysteresis 
curves of phonon gaps in Fig. 13 are reproducible well 
when the starting value of Kb for sweeps is set in the 
small or large Kb regime. Furthermore, even when the 
sweep of Kb is started from the intermediate regime, we 
have observed a tendency that phonon gaps change stay- 
ing on one of the curves shown in Fig. 13 corresponding 
to the increase or decrease in the value of Kb- 



B. Sliding states 

We here investigate the kinetic frictional force of the 
sliding state started from the pinned state. Figure 14 
shows the kinetic frictional force as a function of the slid- 
ing velocity, where the interchain interaction is chosen as 
Kj = 1 . The line with circles shows the result in the case 
that Kb = 2 and K s — 1, which correspond to the valley 
of the maximum static frictional force in Fig. 3(a). For 
comparison, the result for a large Kb limit [Kb = oo), 
where the present model corresponds to the FK model, 
is plotted in the same figure. The frictional force calcu- 
lated by the perturbation theory |||l7j in the case that 
Kb = 2 and K s — 1 is also plotted in this figure. It is 
found that the perturbation theory explains well the nu- 
merical results when the interchain interaction is so weak 
that the maximum static frictional force is vanishing [fl7| . 
For Kb = 2 and K s = 1, the velocity-strengthening and 
velocity-weakening features of the kinetic frictional force 
are recovered and well explained by the perturbation the- 
ory. In the case of the FK model (Kb = oo), however, the 
velocity-strengthening feature in the low-velocity regime 
is destroyed by the large maximum static frictional force 
and the discrepancy between the simulation and the per- 
turbation theory is obvious. Therefore, it is considered 
that the effect of the anomalous pinning in stationary 
states affects the sliding velocity dependence of the ki- 
netic frictional force especially in a low-velocity regime. 

Note here that no hysteretic behavior is observed in the 
velocity dependence of the kinetic frictional force between 
increasing and decreasing processes of the driving force 
when a periodic atomic configuration is employed as the 
initial state of no external force. 

V. SUMMARY 

We have investigated frictional phenomena of an in- 
commensurate two-chain model. By controlling the elas- 
ticity of chains, we have found anomalous pinning be- 
havior, which is accompanied by a successive phase tran- 
sition. The pinned states show complicated behavior 
against the change in elastic parameters and differ ap- 
parently from the conventional breaking of analyticity 
state due to the Aubry transition of the FK model. Un- 
der certain conditions on elasticity and interchain inter- 
action, we have confirmed that the pinning force and the 
maximum static frictional force are anomalously reduced 
or vanishing. Such anomalous pinning behavior is quite 
sensitive to the elasticity of chains, and it is never ex- 
pected for the FK model because the anomalous pinning 
occurs in a characteristic regime where the elasticity of 
both chains becomes important. 

We have also found that the anomalous pinning ef- 
fect affects the kinetic frictional force significantly. Also, 
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in such a pinned state the maximum static frictional 
force is given in the vanishing velocity limit of the ki- 
netic frictional force. In the case that the pinning force 
is weakened anomalously, the velocity dependence of the 
kinetic frictional force shows the velocity-strengthening 
and velocity-weakening features clearly, which are very 
close to those of the kinetic frictional force expected 
in the absence of the breaking of analyticity state. In 
the present paper the overdamped sliding dynamics has 
been investigated in connection with the kinetic frictional 
force. For underdamped dynamics, however, the kinetic 
frictional force may show much more complicated veloc- 
ity dependence, compared with the overdamped case [fl8| . 
It will be discussed in other reports Jl^,|l9| . 

The detailed study of the overall phase diagram in the 
Kb — Kj plane will be reported elsewhere. It might be 
interesting to investigate three-dimensional systems with 
relaxational interfaces. The breaking of analyticity is not 
restricted in one-dimensional systems and is also possi- 
ble in higher dimensions as noticed in Refs. flH- In 
practice, the relationship between the pinning effect and 
the elasticity of lattices in higher dimensions are an in- 
triguing subject in tribology. We may expect more com- 
plicated pinning behavior than that for the present one- 
dimensional model p3). 
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Fig. 1. 

Squared lowest phonon frequency for the FK model 
(K a = Ki = \ and K b = K s = oo). 

Fig. 2. 

Squared lowest phonon frequency vs if;,, (a) Kj = 1 and 
if s = 1, (b) if/ = 0.45 and K s = 1, and (c) if/ = 1 and 
if s = 10. if a is fixed at unity. We will refer to insetted 
numbered arrows later in Figs. 5, 6 and 8. 

Fig. 3. 

Maximum static frictional force vs Kb- (a) Kj = 1 and 
K s = 1, (b) iff = 0.45 and K s = 1. K a is fixed at unity. 

Fig. 4 

Hull functions for if/ = 1, K a = 1, and if 5 = 1. The 
graphs in the left (right) row are the hull functions for 
the upper (lower) chain. The values of Kb, (I) 0.1, (II) 

0. 702, (III) 1.3, (IV) 1.83, (V) 2.11, (VI) 4.09, and (VII) 
8.14, correspond to arrows indicated in Fig. 2(a). 

Fig. 5 

Hull functions for Kj = 1, K a = 1, and K s = 1. The 
values of Kb are (a) 1.18, (b) 1.3, and (c) 1.58, where 
(a) and (c) correspond to local minima of Wj 2 pf and (b) 
corresponds to a local maximum of Wj 2 pf indicated by the 
arrow (III) in Fig. 2(a). 

Fig. 6 

Hull functions for if/ = 0.45, K a = 1, and K s = 1. The 
values of K b , (I) 0.01, (II) 0.15, (III) 0.27, (IV) 0.305, 
(V) 0.43, (VI) 3.21, and (VII) 9.71, correspond to arrows 
indicated in Fig. 2(b). 

Fig. 7 

Hull functions for if/ = 0.45, K a = 1, and K s = 1. The 
values of K b are (a) 0.8, (b) 0.9, and (c) 1.0. 

Fig. 8 

Hull functions for if/ = 1, K a = 1 and K s = 10. The 
graphs in the left (right) row are the hull functions for 
the upper (lower) chain. The values of Kb, (I) 0.01, (II) 

1, (III) 10, and (IV) 100, correspond to arrows indicated 
in Fig. 2(c). 

Fig. 9 

Local lattice structures of the upper and lower chains for 
Kj = 1, K a = 1, and K s = 1. Local atomic displace- 
ments Sui and 5vj are plotted against periodic regular 
sites ic a and jo, in (a) and (b), respectively. The values 
of K b are (1) 0.702, (2) 1.3, (3) 1.83, (4) 2.11, and (5) 
4.09, which correspond to arrows (II), (III), (IV), (V), 
and (VI) in Fig. 2(a), respectively. 



Fig. 10 

Energy quantities vs K b for K / = 1 , K a = 1 , and K s = 1; 

(a) S l a _ a , (b) £ i a _b, (c) #cla- S , and (d) .Bint. 

Fig. 11 

Energy quantities vs Kb for if/ = 0.45, if a = 1, and 
K s = 1; (a) £ ola _ a , (b) E eUi -b, (c) £ c ia- s , and (d) E int . 

Fig. 12 

Energy quantities vs if, for if/ = 1, if a = 1, and 
K s = 10; (a) Bda-a, (b) £ i a _ b , ( c ) ^eia-s, and (d) E int . 

Fig. 13 

Squared lowest phonon frequency vs Kb for if/ = if s = 1 . 
The solid lines with circles indicate the results obtained 
for the sweep with (a) increasing and (b) decreasing val- 
ues of iffc. For comparison, the results in Fig. 2(a) are 
also plotted by the dashed lines with squares. 

Fig. 14 

Kinetic frictional force vs sliding velocity. Triangles de- 
note numerical results obtained in the case that if, — > oo, 
if a = 1, and if/ = 1, i.e., in a case of the FK model. 
Circles and dotted line denote a numerical result and an 
analytic result calculated with perturbation theory, re- 
spectively, for if/ = 1, if a — 1, Kb = 2, and K s = 1. 
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